Domain-averaged Fermi hole versus regional reduced density 
matrices: a critical comparison 
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■ In their recent work Cooper and Ponec [Phys. Chem. Chem. Phys., 2008, 10, 1319- 

S ■ 1329] proposed a "one-electron approximation to domain- averaged Fermi hole (DAFH)" 

' used in electronic population studies. The goal of this comment is to note that the pro- 

posal had been already published within the framework of domain-restricted reduced 
lO ! density matrices (fi-RDM) and to show that it cannot conceptually be considered as 

"^ ' an approximation to DAFH as the authors invoke. 
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The domain- averaged Fermi hole (DAFH) approach to analyze electronic structures has 
been widely used in the recent past^. However, its correlated counterpart (correlated 
hole) has not been considered on an equal footing due to the high computational cost 
required for the use of second-order reduced density matrices (2-RDM), ^D. In ref. 1, 
the authors mention their early work using these matrices^ but omit to advise that 
the first attempt to extend domain-averaged Fermi hole analysis to the correlated case 
has been published by some of us in ref. 3, making explicit use of the ^D structure 
in matrix form. The aim of the authors of ref. 1 is to approach the DAFH analysis 
by a simple model that avoids the use of ^D. The goal of this report is twofold. On 
the one side to show that this one-electron model the authors claim for themselves in 
this reference is nothing but a model previously developed by us within the scenario of 
the domain-restricted first-order reduced density matrices'^, ^D(r2). On the other side 
we attempt to clarify the physical and mathematical reasons, misunderstood in ref. 1, 
showing that what is there called one-electron approximation to DAFH is unsuitable 
for such a task. 

Let us first show that the approach called pseudo-DKFYi [cf. eqn (2.9) of ref. 1] 
is nothing but the model we called symmetric approach to ^D(r2) in ref. 4. Following 
the authors in ref. 1, we will focus attention on closed-shell wave functions. Eqn (11) 
in ref. 4 reads, 

'D}(^) = E ('D^)l. S(fi)f(^D^);. (1) 

k,l 

where (^Dz)* mean the elements of the positive square root matrix arising from the 
spin-free first-order reduced density matrix (1-RDM), ^D, and 8(^2)* =< i\j >q (a* 
in ref. 1) are the overlap integrals over Bader domains Q in an orthonormal basis set 
labeled as i, j, k, I, .... ^. This ^D(r2) arises from a partitioning of ^D in terms of Bader's 
regions, so that ^D* = ^j^ ^DUQ). The expression of eqn (1) in the natural basis set 
in which ^D* = riiSl where Ui and 6^ stand for the natural occupation numbers and 
the Kronecker delta, respectively, yields the symmetric form G^ = nz S{Q) n2 of 



ref. 1 [cf. eqn (2.9)]. Let us note that eqn (1) formulates a true domain-restricted 
first-order reduced density matrix associated with the domain Q, i.e., a representable 
matrix, ^D(r2)^, which is vahd for any type of wave function^, such as independent 
particle models and correlated particle models with non-integer occupation numbers^. 
Our second aim is to show the inconsistencies of the proposal of one-electron 
approximation to DAFH in ref. 1. The errors arise from wrong assumptions of mathe- 
matical and physical nature performed to state the there called pseudo-DAFli^ model 
as an approximation to DAFH. Let us first explain the scenario of the domain-restricted 
decompositions of the 1-RDM, in order to analyze these assumptions. A reduced den- 
sity matrix must fulfil necessary and sufficient conditions to be representable, i.e., to 
assure that there exists a wave function or a statistical ensemble, though unknown, from 
which it derives. The 1-RDM is hermitian, positive semidefinite and bounded^ and is 
normalized to the number of electrons, N. The necessary and sufficient conditions to 
be fulfilled by a closed-shell system domain-restricted first-order density matrix ^D(r2) 
(for any domain Q in the system) are that its eigenvalues, nf must lie within the real 
interval [0,2], i.e., <nf <2^; note that the domains Q define open systems which 
necessarily become described by grand-canonical ensembles^. Let us now consider the 
one-electron approximation to DAFH in ref. 1. This approximation is nothing but 
our symmetric approach to ^D(r2) as we have shown above. Consequently, it turns out 
to be a true domain-restricted first-order density matrix, according to the enumerated 
represent ability conditions. An appropriate model must predict approximate values of 
a determined quantity keeping its physical nature; that is not the case of the symmetric 
model of ^D(i7) in relation with the approximation of DAFH matrices. Thus, as has 
been shown, both locaF and non-local (or integrated) matrix formulations^ of the corre- 
lated hole (DAFH) are not positive semi-definite^'^. Consequently, negative eigenvalues 
(populations) can appear from both formulations, which make this quantity physically 
unacceptable as a density. Furthermore, populations greater than 2 do not fulfilling 
the Pauli principle arise. Examples showing populations out of the interval [0, 2] may 
be found in both the local and domain-averaged formulations of correlated hole^'^. 



Simple systems like N2 molecule at equilibrium geometry calculated at a single-double 
configuration interaction level with the PSI 3.2 package* in a 6-31G basis set, exhibit 
negative eigenvalues (~ —0.03)^. This is a significant negative value which must be 
neglected in DAFH analysis, as was made in refs. 1-3 and consequently the isopycnic 
transformation^" used to localize the eigenvectors of DAHF in those references is no 
longer valid. Also, to neglect the negative eigenvalues of the DAFH matrices permits 
the density to be delocalized due to the particle conservation, i.e., the density integra- 
tion over a domain does not keep the right population for the domain and thus for the 
whole system^. This is not the case of the symmetric model of ^D(r2) which do not 
behave in this way, since it can properly support an isopycnic transformation because 
all its eigenvalues are positive. Thus, this difficulty is completely avoided within the 
framework of domain-restricted first-order reduced density matrix theory, which pro- 
vides a more localized picture of electron distributions as shown in the comparison of 
both ^D(r2) and DAHF pictures^. These results may not be considered as unexpected 
because, as shown above, correlated DAFH are not true particle densities. Namely, 
DAFH may not satisfy the rigorous conditions to be a true first-order reduced density 
matrix, and thus it may not properly describe an open system. Consequently, DAFH 
are not density matrices but different entities and it is a severe conceptual shortcoming 
to approximate DAFH by means of any approach to ^D(r2) including the symmetric 
one. In fact, the above commented results induced us to study decomposition schemes 
in the physical space providing the ^D(r2) matrices^'^. 

It may be noted that DAHF model turns out to be a true density matrix only for 
closed-shell wave functions having 2 or orbital populations. In this case it is equivalent 
to our symmetric approach because then the density cumulant terms of ^D vanish and 
^D is duodempotent (^D = 2 ^D)^. Other important remark to be made in relation with 
the topics treated in ref. 1 is that this one-electron symmetric approach ^D(r2) to DAFH 
has not a common basis of eigenvectors for all the Q domains in the system, i.e., the 
corresponding matrices are not diagonal in the same basis set^. Therefore, calculating 
quantum chemical descriptors depending on the eigenvalues and eigenvectors of two 



domains, as made in refs. 1, 2 and 4, is neither mathematically nor physically rigorous®. 
This drawback may be avoided using a model with a common eigenbasis to all domains 
such as the isopycnic approach to the domain-restricted decomposition of the first-order 
reduced density matrix®. 
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